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Abstract, fn this paper we make a quantitative study of 
the hypothesis that the diffuse Ha emitted from the discs 
of spiral galaxies owes its origin to the ionizing photons 
escaping from H ii regions. The basis of the models is the 
assumption that a fraction of the Lyc luminosity from the 
OB stars within each H ii region escapes from the region, 
leaking into the diffuse gas. A basic input element of any 
such model is a position and luminosity catalogue in Ha 
of the Hii regions in the galaxy under test, down to a 
low limiting luminosity, and we had previously produced 
a catalogue of this type for NGC 157. An initial family 
of models can then be generated in which the Lyc escap- 
ing from an H ii region is parametrized in terms of the 
observed Ha luminosity of the region and the escaping 
fluxes allowed through the diffuse disc gas. These models 
can then be refined using a measured map of H i surface 
density to effect the down-convertion of the Lyc to Ha. For 
NGC 157 an H i map was available. Although its moder- 
ate angular resolution did limit the accuracy with which 
we could test our models, the predicted diffuse Ha surface 
brightness distributions from our models were compared 
with the observed distributions showing that, in general 
terms, the hypothesis of density bounding for the H ii re- 
gions allows us to predict well the spatial distribution of 
the diffuse ionized gas. In the model yielding the best fit 
to the data, the regions of lower luminosity lose a constant 
fraction of their ionizing flux to their surroundings, while 
for H II region luminosities above a specific transition va- 
lue the ionizing escape fraction is a rising function of the 
Ha luminosity. 
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1. Introduction 



Although the diffuse component of Ha emission from a 
disc galaxy has low surface brightness, the integrated Ha 
luminosity emitted by a galaxy in this component is relati- 
vely high, of order 50% of the total radiated by the galaxy, 
for late-type spirals (e.g. Ferguson et al. 1996; Hoopes et 
al. 1996). In Zurita et al. (2000, hereinafter Paper I) we 
made a detailed study of 6 galaxies in which we quanti- 
fied this diffuse component and its geometrical distribu- 
tion across the face of each galaxy. Supporting an obser- 
vation previously noted by Walterbos & Braun (1994) and 
Ferguson et al. (1996), we found a very clear spatial cor- 
relation between the H ii region emission in Ha and that 
of the diffuse emission, leading to a prima facie case for 
considering that OB associations are the sources or, at 
least, the principal sources of the ionizing photons which 
eventually reach and ionize the diffuse gas in the disc as 
a whole. Quantitatively, it appears that the stars in these 
associations can emit sufficient photons above the Lyman 
limit to maintain the diffuse component ionized, as well 
as the gas within the H il regions of the galaxies observed, 
and are an obvious choice for the primary cause of the 
diffuse Ha. They are clearly collectively more luminous in 
the Lyman continuum (Lyc) that the population of white 
dwarfs, for example. However, it has never been entirely 
clear that a major fraction of Lyc photons can escape from 
the immediate surroundings of an OB association: its H II 
region, and even less clear that an escaping photon can be 
transmitted across distances of order a kiloparsec, from 
the nearest OB associations to zones in the diffuse disc 
gas from which significant surface brightness in Ha is ob- 
served. This difficulty was a key motive in leading Sciama 
(1990) to propose his decaying neutrino model for the ion- 
ization of the diffuse component, although recent EUV 
observations of the interstellar medium (ISM) close to the 
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Sun (Bowyer et al. 1999) appear to have ruled out this 
model. 

In the case of the Galaxy, several models have been 
proposed in which transport of ionizing photons over dis- 
tances of order 1 kpc is shown to be viable, at least in 
the neighbourhood of the Sun (Miller & Cox 1993; Dove 
& ShuU 1994). The fact that locally wc can quantify the 
OB stars, and can map the ionized and neutral hydrogen 
in detail, means that models of this sort are useful since 
they can be well compared with data. 

There are two main physical reasons why ionizing pho- 
tons can escape from Hii regions and transit large dis- 
tances in spite of the largo ionization cross section of neu- 
tral H. The principal reason is the inhomogeneity of the 
ISM. It is well established, for example, that although the 
r.m.s. electron density in a luminous Hii region is of or- 
der 1-5 cm~^ (e.g. Kennicutt 1984; Rozas et al. 1996b, 
1999), the electron density determined locally via emis- 
sion line ratios is of order 100 cm~^ (e.g. Zaritsky et al. 
1994). This is interpreted (Osterbrock & Flather 1959; 
Osterbrock 1989) as reflecting the presence of dense knots 
(or clumps of gas) which arc the sites of virtually all the 
optical line emission from the Hii region, within an am- 
bient volume of much lower density gas. The ratio of the 
volume of dense gas to the total volume of the region is 
embodied in the "filling factor", which may be defined 
simply as the ratio of the two volumes. This rather simple 
scheme of dense clumps in a sparse medium (i.e. the inho- 
mogeneous structure of the ISM) is backed by the general 
theory of thermal equilibrium in the interstellar medium, 
treated classically by Field, Goldsmith & Habing (1969): 
a two-phase model, with a warm substrate at T~ 10000 K, 
and cool dense clumps at T<100 K, then by Cox & Smith 
(1974): a three phase model including also the effects of 
supernovae and a component at T> 100000 K, and McKee 
& Ostriker (1977), who produced a more complete three- 
phase model against which much observational material 
within the Galaxy has subsequently been compared. It can 
be demonstrated (see e.g. Spitzer 1978) that in the ISM 
distinct phases should co-exist in pressure equilibrium, 
while states of intermediate density and temperature are 
unstable and decay rapidly. The scenario of dense clouds 
with Ne ^100 cm"'', surrounded by a much lower density 
substrate, with < 1 cm""^, and infiltrated with high 
temperature fully ionized plasma having < 10~^ cm~^ 
has been well verified by many detailed measurements in 
the local interstellar medium (LISM) within a few hundred 
parsec of the sun (Trapero et al. 1995, 1996). The high 
density clouds contribute a major fraction of the mass, 
but occupy a minor fraction of the volume (Trapero et al. 
1996). Studies by Berkhuijsen (1999) in the Galactic ISM 
show that this type of inhomogeneous structure prevails 
over a wide range of mean ISM density and temperature. 
There is evidence for this structure also in external galax- 
ies, as shown by Braun (1997) for a sample of the 11 near- 
est spiral galaxies beyond the Local Group. The fact that 



an inhomogeneous ISM offers a lower effective volume ab- 
sorption coefficient for Lye photons than a homogeneous 
medium is easy to demonstrate, and a very general result, 
common to many types of transfer problems. 

The second effect permitting a fraction of Lyc photons 
to have free paths on kpc scales is a secondary conse- 
quence of the inhomogeneity of the LISM. Since the re- 
combination rate of ionized hydrogen is proportional to 
the product of the electron and proton densities, i.e. effec- 
tively to or to Np, a given luminosity in Lyc photons 
can maintain ionized a volume of plasma inversely propor- 
tional to the square of the density. To see how this affects 
the range of a Lyc photon, we can express this relation 
in terms of the Stromgren radius (the radius of a fully 
ionized plasma sphere maintained by a specific ionizing 
source). The Stromgren radius must vary as ATe"'/^ (or 
Np '^^^). In an inhomogeneous medium, this radius can be 
used as an upper limit to the mean free photon path, if Ng 
is taken as the electron density of the low density phase of 
the ISM. An 03 star yields a Stromgren radius of 140 pc 
in a homogeneous medium of mean density Np=l cm""^ 
(calculation scaled from values in Vacca et al. 1996). The 
corresponding radius in a medium of A'p=0.1 cm~^ will 
be 650 pc, so that a source equivalent to 4 03 stars will 
yield a radius of ^1 kpc in a medium with this density. 
The net result is that in an inhomogeneous medium an 
ionizing source with a strong fiux effectively "burns" its 
own long free path through the low density phase; this is 
the same physical effect as that analyzed by Miller & Cox 
(1993), though it is certainly applicable to more circum- 
stances than those caused only by supernovae (see Mckee 
& Ostricker 1997) in yielding the initial low density ion- 
ized phase. 

We are aware that there are almost certainly other 
sources of energy input into the DIG, which were cited 
in Paper I. These have been inferred from measure- 
ments of emission line ratios, notably [S ii]/Ha,[N iiJ/Ha, 
[Sll]/[Nll] and [OiiiJ/Ho!, as discussed by a number of 
authors, including Reynolds (1984), Reynolds, Haffner & 
Tufte (1999), Greenawalt, Walterbos & Hoopes (1997), 
Hoopes & Walterbos (2000), Collins & Rand (2000), and 
also from the presence of non-thermal velocity widths 
of emission lines (e.g. by Minter & Balser 1999; Tufte, 
Reynolds & Haffner 1999). In the present context, how- 
ever, we are concentrating our attention on what we be- 
lieve to be the principal ionizing source of the DIG, i.e. Lyc 
escape from H ii regions, and our results tend to confirm 
that this approach is useful in practice. 

In modelling an external galaxy to see whether the 
OB stellar sources can account for the observed surface 
brightness morphology of the diffuse Ha-emitting gas, it 
is not possible observationally to resolve structure within 
the ionized gas (nor in the neutral gas). This means that a 
microscopically realistic model of how the inhomogeneity 
affects the transfer of ionizing radiation cannot be tested 
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fairly. In the present paper we have adopted a macrosco- 
pic approach, assuming large-scale homogeneity, so that 
an absorption coefficient with an eff'ective mean value can 
be applied on scales larger than tens of parsec. The range 
of trial values for this coefficient is based on simplified gen- 
eral transfer calculations with realistic photon fluxes and 
gas densities. This method has enabled us to make useful 
comparisons between theory and observation, which go as 
far as is warranted at present by the data, notably in the 
absence of H i surface density maps at high spatial reso- 
lution. The basic aim here is to go at least far enough to 
show whether or not the hypothesis of the OB associations 
as the chief ionizing sources for the DIG can be discarded. 

2. The data used for the models of NGC 157 

In Paper I we published a careful analysis of the observed 
diffuse Ha from the discs of 6 galaxies, but we have chosen 
here only one of these, NGC 157, to model the DIG surface 
brightness distribution in Ha. This is because it is the 
only one of our objects for which we have an H I map with 
anything approaching fully adequate angular resolution. 

2.1. The Ha image 

The observations for this image, and their analysis, were 
described in Rozas et al. (1996a). The galaxy was observed 
at the 4.2 m William Herschel Telescope on La Palma, 
under good seeing conditions. The seeing in the reduced 
final image was 0.8". A luminosity function in Ha for the 
H II regions was extracted and presented in Rozas et al. 
(1996a), after cataloguing 708 Hii regions in total. The 
criterion for considering an image feature as an H ll region 
was that it must contain at least nine contiguous pixels 
each with an intensity of at least three times the r.m.s. 
noise of the local background (see Rozas et al. 1996a for 
further details). 

In Paper I we showed the surface brightness distri- 
bution of the diffuse Ha component, and inferred its inte- 
grated luminosity, its radial distribution within the galaxy 
and the fraction of the total galactic Ha in diffuse emis- 
sion, as a function of radius and in the galaxy as a whole. 

The method employed to separate the Ha flux emit- 
ted from the H ii regions from that of the surroundings 
has benefitted from our previous cataloguing of the H ii 
regions. In the first place, we constructed masked images 
in which the Ha surface brightness of those pixels occu- 
pied by the catalogued Hii regions was set to zero (i.e. 
we assumed that there is no diffuse emission over each 
Hii region), considering a circular approximation for the 
shape of each H ii region. After this, to take into account 
that H II regions are not truly spherical, a surface bright- 
ness cut-off equal to 73 pc cm~^ was employed over the 
masked images to avoid including in the DIG any resid- 
ual emission from the H ii regions. A lower estimate of the 
total Ha emitted from the DIG of the galaxy was then 



obtained by integration of this masked image. The upper 
estimate was obtained in a similar way, but including an 
estimate for the diffuse emission above each Hii region. 
The detection thresold for the DIG of NGC 157 was 1.0 
pc cm~^. Detailed explanations of the isolation procedure 
and uncertainties in the measurements of the total Ha 
emission from the DIG are given in Paper I. 

As stated in paper I, we have not made any correction, 
either for the dust extinction within the emitting galaxy, 
or for Galactic extinction. 

2.2. The Hi map 

NGC 157 was observed by Ryder et al. (1998), with the 
VLA using the C and D configurations. The on-source 
exposure time with the C array was 3.7 hours, and with 
the D-array 2.8 hours. The two data sets were processed 
separatedly and combined. The total velocity range cov- 
ered was 581 km with a channel increment of 5.2 km 
s~^. Two distinct data cubes were produced, using differ- 
ent weightings in the imaging process (see Ryder et al. 
1998 for a detailed description of the observations, the 
reductions and the processing tasks). Here we used the 
H I surface density map obtained from the cube produced 
using uniform weighting, which yields better angular re- 
solution than natural weighting. The beam width at half 
power for the map used here is 18"xl2" and the r.m.s. 
noise 1.3 mJy beam~^. 

3. Descriptions of the models for the propagation 
of the ionizing photons 

In the present section we explain the simplified theoreti- 
cal models which predict the ionizing photon distribution 
in the disc of the galaxy NGC 157, assuming that they 
are produced within and escape from the catalogued H ii 
regions. Different laws governing the escaping photon frac- 
tion and the absorption as they transit the ISM yield a set 
of models. In these models, each measured Hii region is 
taken as an ionizing source and fraction of the measured 
ionizing flux within each region is taken as the luminosity 
of the escaping flux, available to maintain the ionization 
of the DIG. The Ha surface brightness of the DIG is com- 
puted by allowing the Lyc flux from a source to radiate 
within a thin disc in an inverse square regime of dilu- 
tion. In the simplest case, the fields from all the observed 
H II regions are co-added and the result compared to the 
observed surface brightness distribution. Implicit in this 
procedure is that the DIG surface brightness in Ha is lin- 
early proportional to the underlying Lyc field, which is 
the same as assuming that the disc consists of a slab of 
HI, whose properties do not vary from place to place. It is 
clear that the observed Ha surface brightness must be the 
projection of a three-dimensional DIG column onto the 
plane of the galaxy observed (if it is face-on) . Given that 
the scale-height of the DIG, as shown within our Galaxy 
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(Reynolds 1989) is of order 1 kpc, and that even the largest 
H II regions do not project more than some three hundred 
pc out of the mid-plane, the simplified geometry employed 
here will yield a useful result for the integrated columns 
perpendicular to the plane viewed in a generally face-on 
direction. 

A second step in bringing the model into better rela- 
tion with the observations is to use the observed Hi sur- 
face density distribution to modulate the Lyc flux, taking 
a uniform mass absorption coefficient to represent atte- 
nuation by gas plus dust extinction, and a third step is to 
take the product of the previously computed Lyc field and 
the Hi surface density for the predicted Ha distribution, 
assuming that the resulting surface brightness is the result 
of the irradiation of local neutral hydrogen by an incident 
ionizing field. The results from this type of modelling are 
expected to be useful approximations provided that the 
ionized column density is not a large fraction of the neu- 
tral column density, an assumption which we tested for 
NGC 157, and which holds everywhere in the disc outside 
the Hii regions (see Sec ||). 

An interesting invariance relation which helps to a- 
ccount for the relative lack of apparent structure in the 
DIG, in spite of the assumed inhomogeneity of the clouds 
which compose it, and which also helps to account for the 
fact that an essentially two-dimensional approximation 
can satisfy a three dimensional flux distribution, is that 
the emission measure in Ha emitted by an optically thick 
cloud irradiated from outside is independent of the cloud 
density, and depends only on the intensity of the irradiat- 
ing field. We can show this quite directly. In an optically 
thick case, where all the incident photons are absorbed, 
the intensity of the recombination radiation emitted (in 
the present case Ha) can be related to the incident ionizing 
radiation, using the "on the spot" approximation (Ostcr- 
brock 1989), in which any ionizing photons emitted by the 
ionized gas are reabsorbed locally. In this case, the number 
of ionizations per unit volume of the cloud in which the 
external flux has penetrated will be equal to the number 
of recombinations to excited levels of the hydrogen: 

/•oo 

Nho / F,a,dv^N,NpaB{H\T) (1) 

where the atomic hydrogen density, the electron density 
and the proton density are denoted respectively by Nh° , 
Ne and Np, where aB{H°,T) is the recombination coef- 
ficient to excited levels of hydrogen, F^, is the number of 
incident photons per unit area, time and frequency inter- 
val, and fly is the ionization cross-section for hydrogen by 
photons with energy hv above the Lyman limit hvo ■ If we 
integrate Eq. ^ down to the depth within the cloud within 
which all the incident radiation is absorbed (the optically 
thick case) using the equation of radiative transfer in the 
form: 

dF^/ds = -NHoa^F^ (2) 



we have: 

/ {/ ~dF^} dv = aB{H° ,T) / N^Npds (3) 

Ji^o Jf{s=o) Jo 

The integrated intensity of the emitted Ha from a co- 
lumn of length L along which the electron and proton den- 
sities are given by N^. and Np respectively is proportional 
to the integral on the right hand side of Eq. || conventio- 
nally termed the emission measure, EM, of the column. 
We can rewrite Eq. |^, incorporating the appropriate nu- 
merical constants, in the simpler form; 

F{s = 0) = .4 aB{H°,T) EM(in pc cm"^) (4) 

where F{s — 0) is the incident flux in the surface of the 
cloud, A is a constant, whose numerical value for F{s = 0) 
in units of phot, s'^cm'^ is 3.08 x 10^^ and aB{H° ,T) 
does not show strong dependence on T in the range nor- 
mally found in the warm ionized ISM (as ~ 2.59 x 10~^^ 
cm'^ s~^). 

Eq. ^ implies that irrespective of the total column 
depth of the cloud, if it is optically thick to ionizing ra- 
diation it will develop an emitting skin for which the Ha 
intensity will be proportional to the ionizing flux. Thus, 
seen from above, a set of H i clouds with varying hydrogen 
column density bathed from the side in uniform ionizing 
radiation will all show the same surface emissivity. This 
result accounts in large part for the uniform appearance 
of the diffuse Ha from a galactic disc, in spite of the in- 
trinsic dumpiness of the substrate. This result does not 
hold for an optically thin cloud, so that the net diffuse 
emission from a projected disc should not be entirely free 
from local variations. The result above also implies that, 
if the ionization of the DIG is indeed due to OB stellar 
clusters, the observed mean EM in the disc relatively far 
from ionizing sources should show a rough proportionality 
to the Hi column density. This is because if the Hi is in 
discrete clumps, optically thick to Lyc (but not in general 
to Ha) , the integrated EM averaged over a column should 
be proportional to the integrated surface area of clumps 
per projected column cross-section, which will be a rising 
function of the integrated column density in HI, although 
the exact relation will depend on the details of the clump 
distribution. These considerations for the relation of the 
observed EM of the DIG to the gas distribution in the 
disc are valid for whichever detailed hypothesis we adopt 
about the relation of the observed Ha in the H ii regions 
to their ionizing luminosities in Lyc. 

3.1. Constant escape fraction (case A) 

The simplest assumption to make when considering Hll 
regions as the zones of origin for the photons which io- 
nize the DIG is that a constant fraction of the Lyc from 
the OB stars in any Hii region escapes from the region 
into the surrounding diffuse gas. A number of previous 
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studies have suggested that a significant fraction of the 
total ionizing luminosity of OB stars can and does escape 
from their surrounding H il regions, and is available to io- 
nize the diffuse ISM. This conclusion is supported by the 
geometrical correlation observed between the sites of H il 
regions and the surface brightness of the DIG (Ferguson 
et al. 1996). The escaping flux from a given H ii region was 
not quantified by Ferguson et al. (1996) who did not pro- 
duce a photometric catalogue of Hii region luminosities, 
though their estimate of the total flux required to ioni- 
ze the DIG was comparable to our estimates in Paper I, 
i.e. of order 50% of the total for the whole disc. McKee 
& Williams (1997), estimated that approximately 30% of 
the ionizing photons observed within the H ii regions of 
the Galaxy, defined by the extent of their cm wavelength 
continuum emission, in fact escape from the Hii regions. 
In general, the number of ionizing photons absorbed in an 
Hii region observed in Ha will be higher than that ab- 
sorbed only within the bright central zones, which define 
the region at cm wavelengths, but it is not apparent how 
to define reliably a core-halo boundary in Ha. We chose a 
30% constant escape fraction as our starting point for the 
"constant escape fraction models" of NGC 157. In these 
models, each measured Hli region is taken as an ionizing 
source and a constant fraction of the measured ionizing 
flux within each region is taken as the luminosity of the 
escaping flux available to maintain the ionization of the 
DIG. The Ha surface brightness of the DIG is computed 
by allowing the Lyc flux from a source to radiate within 
a thin disc in an inverse square regime of dilution. 

3.2. Density bounding at high luminosity (case B) 

In Paper I we outlined a possible model for the escape of 
ionizing photons from Hii regions, in which essentially all 
the Lyc escapes from the regions of highest luminosity, i.e. 
those with observed luminosity in Ha, hna, greater than 
a critical transition value l^str (the "Stromgren transi- 
tion luminosity"). This follows from a series of studies in 
which evidence has accumulated pointing to the fact that 
H II regions of very high luminosity are density bounded, 
and that an increasing fraction of their Lyc fluxes escape 
as their luminosities increase (Rozas et al. 1998; Beck- 
man et al. 2000) . We have developed a general parametric 
formulation showing the dependence of the escaping ion- 
izing luminosity as a function of the physical parameters 
of the Hii regions with Ha luminosity greater than l^str, 
whereas no escape of Lyc is assumed for Hii regions in 
which IjHa <L5tr in this case B. 

Let N/za be the number of ionizing photons per unit 
time from OB associations which ionize the surrounding 
cloud and are down-converted in Ha emission within an 
H II regiorj^. Assuming that an H li region contains gas in 

^ The relation between Ni/^ (number of ionizing photons 
per unit time from the OB associations down-converted in Ha 



small condensations or clumps with electron density pc, 
^Ha can be expressed as: 

- ynalVe^</.R3 (5) 

where a^^ is the effective Ha recombination coefficient, 
(j) is the filling factor (or fraction of the total volume oc- 
cupied by the condensations, as explained above), n is the 
number of Lyc photons necessary to yield one Ha photon 
(n=2.233 Lyc photons/Ha photon for T~10'' K; Oster- 
brock 1989) and Rs is the Stromgren radius of the Hii 
region. 

The ionized mass of the Hii region can be written as: 
= f mH4>PcdV = — ^ (6) 

where mu is the mass of an atom of hydrogen and the 
integral is performed over the entire volume {V) of the 
Hii region. 

Density bounded H ii regions are fully ionized, so for 
these regions the mass of the placental cloud (Mcz) is equal 
to the mass of ionized hydrogen in the H ii region. Com- 
bining Eqs. 1^ and ^ and normalizing the quantities (which 
will be noted with " ) by their value at the Stromgren tran- 
sition, we find that 



Nff„ = -^-^^ = , 7 , " = Me/ (7) 

\^Ha)StT (Pc)str V^cljStr 

This equation holds for density bounded Hii regions. 

In Beckman et al. (2000) a simplified scaling calcula- 
tion is presented, in which the required relation between 
cloud-mass and the ionizing luminosity (L^) emitted by 
a young cluster of stars which has formed in a placental 
cloud of mass Md is given by 

L, = A;(M,)" = fcr(Me/)"% (8) 

where fc, j are constants that depend principally on the 
metallicity (Beckman et al. 2000) and M^ is the total ste- 
llar mass in the OB association. The index a represents 
a parametrized relation between Lj and the stellar mass 
{M^,) in the cloud, while the index e is the parametrized 
relation between the the stellar mass and Md . The condi- 
tion that the rate of production of ionizing photons fully 
ionizes and overflows the placental cloud implies that ae 
should be greater than unity (see Beckman et al. 2000 for 
details). Observational evidence showing that ae should 
indeed be greater than unity was reviewed in Beckman et 
al. (2000), and a value close to 1.6 was adopted as the 

emission) and the observed luminosity in Ha of an Hll region, 
Lhu is N_f/ci = nLHa/(hvHa), where hvHa is the energy of an 
Ha photon and n=2.233 Lyc photons/ Ha photon for T^IO* 
K (Osterbrock 1989). 
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Fig. 1. Fraction of the total ionizing photon flux within an 
H II region, which escapes from it in the mixed model C 
(see text for details). The relation of fcsc (in the figure) to preceding Sec. |31 
/ (Eq. |l^) is given by fcsc=//(l+/). Approximately 80% of 
the total ionizing radiation escapes from the most luminous 
H II regions in NGC 157 in this model (i.e. nearly 4 times the 
ionizing luminosity down-converted to Ha emission within 
the region). 



from the H ii region. Note that if / (Eq. |T^) is the fraction 
of ionizing photons down-converted to Ha (Nho) which 
escapes from the Hii region, the relation between the two 
quantities is given by fosc=//(l+/)- 

In the models, to maintain the family of calculated 
cases as simple as possible, we maintained ae at a con- 
stant value of 1.6, and pc constant, varying in order to 
reproduce the range of observed surface brightness in Ha 
V. Lffa found for the regions in NGC 157 by Rozas et 
al. (1996b). We note that at the highest luminosities the 
escape fraction rises to values well above 0.5, which means 
that more than 50% of the Lyc photons escape from the 
region. We can use the interpolated curve in Fig. ^ to 
compute an escape fraction from each Hil region in the 
catalogue of regions for the galaxy, so that the absolute 
luminosity in ionizing photons escaping from each region 
is then tabulated, and can be used as input to model the 
production of the diffuse Ha in the same way as for the 
models with constant escape fraction, as explained in the 



3.3. Mixed models (case C) 



In this model we assume that the major contribution of 
ionizing photons available to ionize the DIG comes from 
the most luminous H ii regions (those with log l^Ha ^38.6, 



best estimate. Combining Eqs. 
to (Li)str: 



and 1^ and normalizing 



(9) 



Under these conditions we can compute the escape frac- 
tion (/) of ionizing photons down-converted in Ha iJ^Ha) 
which escapes from the most luminous regions as a func- 
tion of the observed luminosity in Ha: 



(i^Ha) Pc — 1 



(10) 



To estimate / it is necessary to know how the clump den- 
sity varies with the H II region luminosity. Measurements 
of Ha surface brightness alone will not allow us to detect 
whether pc is increasing as the luminosity increases, but 
do tell us (see Rozas et al. 1996b, 1998) that the prod- 
uct pc4>^/'^, which is equal to < Ne >rms, increases by a 
factor close to 2 between the Stromgren luminosity and 
the luminosities of the brightest regions, and enables us 
to use Eq. |l^ to set bounds to the predicted escaping pho- 
ton flux from density bounded regions, as a function of 
luminosity. An example of this is given in Fig. (for lumi- 
nosities greater than l^str) in which the two limiting cases 
of varying pc and fixed 4> or varying and fixed pc are 
taken and used to predict the escape fractions of ionizing 
photons versus luminosity. Fig. Q shows fcsc , which is the 
fraction of the total ionizing photon flux which escapes 



in erg s~^), whose escape fraction is given by Eq. 10, but a 
minor contribution escapes from the less luminous regions 
(with log IjHa < 38.6, in erg s~^), for which a constant 
escape fraction of ionizing photons is assumed (/=20%). 

Inspection of the Ha image of NGC 157 showed us 
that the diffuse component is relatively strong near an 
H II region, and that this observational correlation appears 
to hold not only for the most luminous regions, but also 
for regions with luminosities below hHa=T-istr- The impli- 
cation is that even those regions with lower luminosities 
are leaky. Such an inference is entirely plausible given the 
clumpy structure of any interstellar cloud and, in particu- 
lar, of the clouds which make up the Hii regions. Proof of 
the clumpy nature of the regions is found by comparing 
the r.m.s. electron density in a region (Rozas et al. 1996b, 
1999, 2000) obtained from a diametral emission measure in 
Ha with in situ measurements of electron density obtained 
via line ratios (Zaritsky et al. 1994), as explained above in 
Sec. |l|. Relatively crude models using this structure give 
as a result an escaping fraction for the Lyc produced by 
their OB populations, a fraction which should vary little 
if the density structure of a region depends little on its 
luminosity. Measurements of filling factors show that, at 
any rate for Lnct < ^str this condition appears to hold 
rather well. However, as explained above in Sec. 3.1 at 
the highest luminosities the r.m.s. electron density shows 
higher values (Rozas et al. 1996b, 1999, 2000) which we in- 
terpret as due to an increased effective filling factor, which 
in turn is explained if the clumps near the stellar sources 
become fully ionized, releasing more Lyc to ionize clumps 
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where {xj,yj) are the coordinates of the region selected. 
We then sum over all the regions to obtain the final ioni- 
zing photon flux density map from the contribution of all 
Hii regions: 



H 11 regions iii Llie deprujecLed disc uf NGC 157. The veilical 
direction of the plot is coincident with the major axis of the 
galaxy. Symbols show ranges of log Lho ■ 



further out within the region. In this model each indivi- 
dual clump behaves as an ionization bounded or a density 
bounded system, according to the ionizing luminosity of 
the region, and the distance of the clump from the star 
cluster. The Hii region as a whole then starts to leak a 
higher fraction of its Lyc luminosity when the inner dense 
clumps become density bounded, which will occur for re- 
gions of high stellar luminosity. Without going into further 
detail, we have taken this into account empirically in our 
termed "mixed models" by using a constant escape frac- 
tion for IjHa < Lstr s-nd an increasing escape fraction for 
Lifa > Lstr, following Fig. 0. 

4. Comparing observations with the model 
predictions 

Models with no computed in-plane extinction 

The basis for all the models is the positional and luminosi- 
ty catalogue of all the detected Hii regions in NGC 157^. 
In Fig. ^ we show this in display form, i.e. the positions 
and luminosities of all the regions used in the models are 
indicated in this map deprojected into the plane of the 
sky. To use the regions as sources for the diffuse compo- 
nent, we first apply the escape factor appropriate to the 
model: constant factor (case A), pure density bounding 
above LHa=^str (case B), or mixed model (case C), to 
each and every one of the regions. Treating each region 
as a point source with effective Lyc luminosity equal to 
that corresponding to its measured Ha output, multiplied 
by the escape factor, we then model the diffuse emission 
using a set of assumptions of increasing complexity. 

In this section we assume a uniform HI layer of con- 
stant thickness. In the simplest case, we compute a flux 
density Ij{x, y) at a point (x, y) in the disc due to a region 
(j) emitting Dj escaping photons per second, via: 



lj{x,y) 



At:[{x-XjY + {y-yjf 



(11) 



N 



Hx,y) = ^/j(x,?;) 



(12) 



^ For the criteria of the lower limit of detectability the reader 
is referred to Rozas et al. 1996a. 



where TV is the total number of catalogued regions (weaker 
regions would make a negligible impact on the result using 
any of the assumptions taken in this work (and can gen- 
uinely be neglected). A sample result of computing /(x, y) 
in this way is shown in Fig. |3^ where we have taken the 
simplest case, i.e. a hypothetically constant fraction of es- 
caping photons for the Dj values. We have compared this 
computation directly, in Fig. ||(right) with the deprojec- 
tion into the plane of the sky of the diffuse Ha measured 
in NGC 157, and we can see that the result is to say the 
least encouraging. All the broad intensity and morpholog- 
ical features of the diffuse Ha are reproduced, although 
inspecting the fine details we can pick out certain dis- 
crepancies. We should state now that in the rest of the 
discussion rather than deproject the observed diffuse light 
distribution into the plane of the sky we will project the 
model into the plane of the galaxy disc, as this is an o- 
peration which of course is not affected by low signal to 
noise ratio. 

The fact that the morphological agreement between 
model and observations is so good is perhaps a little sur- 
prising, since the conditions imposed are so simple. With- 
out further refinement we can use this result to criticize 
the hypothesis that the diffuse Ha is originated by de- 
caying massive neutrinos, although this model by Sciama 
(1990) has been already found wanting via direct EUV 
observations (Bowyer et al. 1999). We can do this by com- 
paring the measured distribution of H i column density in 
Fig. ^with the Ha distribution and model prediction. The 
Hi density contours for NGC 157 (Ryder et al. 1998) have 
been overlaid on the Ha image of the galaxy. The angular 
resolution is poor by Ha standards, but we can see clearly 
that zones of high surface brightness in diffuse Ha do not 
coincide with zones of high Hi column density, a result 
which contradicts the predictions of the massive neutrino 
decay hypothesis. This lack of correlation is maintained if 
we smooth the Ha image to match the resolution of the 
H 1 image. The leaky H ii region model, even in its simplest 
form as presented here, offers a much better solution. 

In Fig. ^ we present the results of applying the sim- 
ple geometrical propagation model to leaky H ii regions in 
three regimes: constant escape fraction, density bounding 
only for L^q > Lstr, and the "mixed model" as descri- 
bed above, compared with the observed distribution of the 
diffuse Ha. It is clear that the model with density bound- 
ing only for Lhu > Lstr- yields a worse fit to the surface 
brightness than the other two. On the basis of this re- 
sult we can already begin to question the hypothesis put 
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Fig. 4. Colour representation of the modelled ionizing photon density in the disc of NGC 157, assuming that this 
distribution is produced by the escaping photons from the most luminous regions alone (case B, a.), from all the 
catalogued Hii regions (case A, b.) or from all Hii regions but with different escape fractions (case C, c). The images 
have been normalized to the mean value in the ellipse determined by the observed Ha emission of the galaxy (see 
Paper I). The observed DIG is shown in d. for comparison (also normalized). 
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RA (2000.0) 

Fig. 5. Ha surface brightness density of the observed DIG 
in NGC 157 with Hi column density contours superposed. 
The density contours are 1.45, 2.90, 4.35, 5.80, 8.70, 11.60 
Mq/pc^ (col umn density increasing from finer to thicker 
contours). The two quantities (Ha surface brightness and 
H I density) should be proportional in the "Sciama neu- 
trino" model. We can see that (in spite of the difference 
in resolution between the two images) the H i distribution 
does not correlate at all with the Ha emission by the DIG. 
Darker areas in the image correspond to higher levels of Ha 
emission by the DIG of the galaxy. Thicker contours repre- 
sent the higher Hi column densities. The parallel shading 
in the corners of the image is an artifact generated in the 
rotation of the image to obtain the correct orientation (the 
same effect is seen in Figs. |7|i, || and 



forward in Beckman et al. (2000), that only those regions 
with Lho! > Lstr are density bounded. This is particularly 
clear for the diffuse emission from the group of H ii regions 
forming the uppermost arm of NGC 157. Although most 
of these regions have Lho < Lstr the arm is surrounded 
by diffuse emission, which is not well reproduced in the 
"density bounding" model where Lyc photons escape only 
from the regions with Lhu > Lstr- In all the models the 
distribution of the diffuse emission is less confined to the 
zones around the arms than in the galaxy as observed, and 
this leads us to conclude that the underlying assumptions 
are rather too simplified. 



4-. 2. Models incorporating the observed Hi column 
density distribution 

One of the simplifying assumptions implicit in the results 
of Fig. ^ is that the H i which down-converts the escap- 
ing Lyc photons to Ha is in a disc of uniform thickness. 
We carried out a check to see whether at least part of 
the discrepancy between models and observations could 
be attributed to this. In Fig. ^,b,c, we show ratios of the 
modelled Ha surface brightness to the observed surface 
brightness, for the three cases presented in Fig. ^. The 
ratios are those of Figs. ^a,b and c, respectively, divided 
by Fig. and smoothed to four times the original pixel 
size. In Fig. ^ we have shown the measured H i surface 
density for comparison. We note that the zones where the 
ratio of predicted to observed Ha is high in Figs. ||a,b,c, 
coincide with the zones where the observed Hi column 
density is low in Fig. which implies that our neglect of 
the H I column density variations in our model predictions 
of Fig. II needs to be corrected. We have performed such 
a correction, again in the simplest way possible, by mul- 
tiplying directly the predicted Lyc field in a model by the 
Hi column density. The results are shown in Fig. ^ where 
we can see that, globally, the distribution of Ha from the 
DIG is better reproduced in general than in Fig. |[ in that 
the diffuse component is more intense closer to the arms, 
and less so, further away. However in one aspect the "H i 
corrected" models reproduce the oberved DIG less well 
than the previous ones, i.e. in the zone around the centre 
of the galaxy, where the predicted intensity is too low in 
the Hi corrected models. One reason for this is that we 
have simply assumed that the Ha emission ought to be lin- 
early proportional to the Hi column density. This would 
be a fair approximation only if the H ii column density 
were small compared to that of H i; if not it would seri- 
ously under-estimate the predicted Ha. We carried out a 
test for this, illustrated in Fig. ||, in which we first esti- 
mated the H II column density, and then compared it with 
that of H I. Since we measure directly a surface bright- 
ness, i.e. an emission measure (in Ha) our estimate of the 
column density cannot be unique, without a further as- 
sumption, which is the effective thickness (h) of the Ha 
emitting layer. Unfortunately the result depends on the 
square root of this thickness, which cannot be measured. 
In Fig. H we have used 900 pc for this thickness, taking the 
"Galactic Reynolds layer" as our guide. Fig. ^ maps the 
resulting Hii column density in the DIG, and in Fig. ||b 
we have divided this into the observed H I column density, 
yielding a map of Hi/Hii. We see from this map that the 
ratio is between 5 and 15 over most of the galaxy disc, 
but that near the centre it drops to low values, so that 
we would expect our H i-corrected models to fall short of 
the observed surface brightness in Ha near the centre, as 
found in Fig. 0. 

One problem in allowing for the H I in this modelling 
process is that the angular resolution in the H i map of 
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Fig. 6. First three panels: ratio of the modelled photon distribution, for the three cases shown in 
Fig. ^,b,c, to the observed DIG emission in Ha. Pixels have been binned to a new pixel size equal 
to 4. Only the lowest Hi density contours have been superposed (of. Fig. ^ and panel d.) There is 
good spatial coincidence of the zones where the models with an assumed layer of H I of constant 
thickness overestimate the emission (high ratios correspond to darker areas), with the "holes" in 
the measured H i distribution. This distribution is shown in grey scale in the fourth panel, in which 
the darkest areas correspond to the highest H i column densities. 



NGC 157 is no better than 12", which is a factor of 15 
lower than we would wish, in order to make adequate use 
of our Ha map. In spite of this we have gone one step fur- 
ther in our modelling, taking into account the absorption 
by the intervening medium of the propagating Lyc field. 
This is done by substituting Eq. |ll| by 

(13) 

= F^ exp {-Tj) 



in which Fj is the flux incident on unit path length at 
{x,y) from a source j, and Tj is the optical depth along 
the path to {x,y), which is the integrated product of the 
effective volume absorption coefficient (fc) and the volume 
density {pm) along the path between the emitting Hii 
region j and the point {x, y). In implementing this model, 
the effective absorption coefficient is taken as an empiri- 
cal parameter, which includes absorption by the H i (the 
cause of the Ha emission) and also dust extinction, with 
an assumed constant dust to gas ratio. 
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Fig. 7. Monochrome represeiftci,<!!f5fi' of the modelled diffuse Ha emi^^diPTfcr the galaxy NGC 157. 
Each figure represents the product of the photon distribution and the H i density at each point of 
the galaxy for the cases B, A and C shown in Fig. ^ All four figures have been normalized to the 
mean value in the ellipse determined by the observed Ha emission of the galaxy and plotted using 
the same grey scale. The different resolutions between the Ha and the Hi data does not allow us 
to extract a definitive conclusion from these maps, but the broad agreement in the overall DIG 
distributions predicted with that observed, strongly favours the hypothesis that photons leaking 
from Hii regions do power the DIG (see Sec. ^for a discussion). 



Using Eq. we can express the decrement, dij /ds of 
ionizing radiation per unit length ds at the point (a;, y) by 



dl^ dFi / \ 



I CXp ( — Tj ) 

ds 



(14) 



-jf exp(^Tj) + i~kps)Fj 



The first term here represents the decrement due to geo- 
metrical dilution, and the second term the component due 
to absorption, both terms being intuitively clear. In order 
to get a reasonable value of k without three-dimensional 
modelling of the clumpy medium, we take an empirical 
approach based on the observation that the diffuse com- 
ponent extends to distances of order 0.5 kpc from the con- 



centrations of Hii regions. With this starting point we take 
a trial value of k as that which would reduce a flux to 1% 
of its initial value over a range of 500 pc. Taking an a- 
verage Hi scale height of 100 pc as a constant value over 
the disc and a mean column density in the disc as that ob- 
served: 8.6 X 10^° cm~^, we find a value for k of 5.3x10"^^ 
cm^, or more generally 5.3xl0~^^ x zioo cm^, where zioo 
is the value of the scale height in units of 100 pc. In fact, 
this value for k ought to be an upper limit, because the 
condition that the flux be reduced by a factor 100 over 
a distance of 500 pc is a stronger condition than that re- 
quired by the empirical map of the DIG. In our numerical 
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Fig. 8. a. H+ column density map assuming an effective thickness of the Ha emitting layer {h) 
equal to 900 pc. Since the ionized mass per unit area is given by: m(H"'") =< Nf, >rms mnh and 
the emission measure EM=< >rms the relation between the number of ionized atoms of H+ 
per cm^, n(H"'"), the EM and h is given by: n(H+) ~ VEM/i^/^. The mean value for the DIG is 
~ 5 X 10^° H+ ions cm~^. For other scale heights the monochrome bar calibration changes by a 
factor ( goope )^^^' R-^-tio of the Hi column density to the H+ column density assuming an H+ 
effective layer thickness equal to 900 pc in the disc of NGC 157. In spite of the different image 
resolutions, in the DIG area most distant from Hii regions (where surface brightness gradients 
are nearly constant, see Paper I) the figure will give a reasonable estimate of the ratio of the two 
column densities. In these areas, the Hi column density is between 5 and 12 times greater than 
the H+ column density. 



models we took values of k ranging between 5.3 x 10 
and 5.3 x 10~^^ x ^loo cm^. 

Our models are of two kinds: mixed models in which 
the absorption coefficient k is allowed to vary, but in which 
the H I is assumed to take the form of a disc of constant 
thickness, and models in which an optimized value for k 
is combined with the observed Hi column density map. 
The two types of models are illustrated in Figs. ^ and ^ 
respectively. In Fig. ^ we can see that the best approxima- 
tion to the observed DIG surface brightness distribution is 
found using an empirical coefficient k = 5.3 x 10~^^ x zioo 
cm~^, and in Fig. |l^ we can see the results of using this 
value with two different assumptions about the Lyc escape 
fraction: a constant escape fraction and a constant fraction 
below Lhu = Lstn and a rising fraction for higher H ii re- 
gion luminosities (cases A and C respectively) . The results 
shown in Figs. ^ and |l^ illustrate the fact that we are find- 
ing it difficult to obtain real improvements in our results 
with increasingly refined assumptions. Perhaps Figs, 
y between them give the closest reproductions of the 



observations (e.g. Fig. ||d), and seem to indicate that a 
constant escape fraction for all the regions gives the best 
fit to the observed DIG. However there are a few necessary 
cautionary statements to make. We can see that Fig. 
in which the fullest set of assumptions has been brought 
in: Lyc field modified by absorption, and emission in Ha 
proportional to the H i column density, gives a clearly in- 
ferior reproduction of the DIG morphology than Fig. 
assuming a constant H i column density. This is due to two 
effects, firstly the H i column density should not be used in 
this way where the Hil column is a major fraction of the 
total H column, as occurs near the centre of NGC 157 (see 
Fig. H and discussion in Sec. ||). Secondly, we are coming 
up against the limits imposed by the relatively poor an- 
gular resolution in the Hi map. This is not so important 
far from the sources, where the Lyc field gradient is low, 
but is critical close to the Hll regions, and means that 
models with any further degree of refinement cannot be 
properly empirically tested. The combination of these two 
effects can be already leading to misleading conclusions in 
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Fig. 9. Same as Fig. Q but for models including absorption by neutral hydrogen. In this case, all figures assume 
case C (mixed model) for the escape fraction of ionizing photons. The effective volume absorption constant for 
each case is shown in the left top corner of each representation. 
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Fig. 10. Same as Fig. but for models including absorption a. for the mixed model and b. for the constant 
escape fraction model. 
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terms of the fine detail. An apparent defect in those mod- 
els where the high luminosity regions have an increasing 
Lye escape fraction is that the predicted Ha near these 
regions is higher than that observed. Such an effect would 
be expected if the H ii/H i column density near the regions 
was raised (the same effect that we see close to the cen- 
tre of the galaxy in the constant Hi thickness model of 
Fig. ^,c) due to the ionization itself. The effect could be 
detected only with Hi resolution of order 1", i.e. compa- 
rable to that in Ha. In the absence of such data we are 
in no position to detect the effect, so the impression that 
the constant escape fraction models are better than the 
mixed models could well be an incorrect interpretation. 

5. Discussion and conclusions 

We have used a modelling technique to examine the hy- 
pothesis that the Hii regions in disc galaxies are suffi- 
ciently leaky that they allow the ionizing photons from 
their OB stellar populations to escape with intensity high 
enough to produce the diffuse Ha observed across the disc 
of a galaxy. We have used NGC 157 as a test object here, 
because one of the observational elements required in the 
models is a map of the H i column density of the object 
to complement an existing map in calibrated Ha surface 
brightness. A second element in the modelling process is 
to derive a catalogued intensity calibrated map of the indi- 
vidual Hii regions in Ha, complemented by the measured 
surface brightness over the full diffuse component. These 
data are not trivial to acquire, and for the galaxies for 
which we had carried out calibrated Ha mapping of the 
necessary precision, only NGC 157 had an available map 
in Hi of anything approaching adequate angular resolu- 
tion. In spite of the approximations and data limitations 
mentioned freely at the relevant points in the text, we have 
found a very fair coincidence between the distribution of 
ionizing photons, and the H+ distribution across the face 
of NGC 157, as shown in Figs. |, § |[ and |l|. This leads 
us to conclude that in spite of the simplicity and the ap- 
proximations with which the models have been generated, 
and the low angular resolution of the H i data, the hypoth- 
esis that the ionization of the DIG in NGC 157 is caused 
by the Lyc photons escaping from the Hii regions does 
receive major support from this exercise. 

The fact that the hypothetical model based exclusively 
on the Hi distribution, and designed as a test of the 
Sciama neutrino decay theory (Fig. ^) gives such com- 
plete disagreement with the observed Ha distribution in 
the DIG not only serves as a further rebuttal of this the- 
ory, but points up the rather good agreement of the models 
using leaky H II regions as the ionizing sources. 

In order to estimate quantitatively the goodness of the 
models we performed a series of numerical tests. These 
consisted of an analysis of the residuals obtained when 
subtracting the observed Ha emission of the DIG from the 
modelled emission (both observations and models normal- 



Table 1. Lyc continuum photon flux escaping from the 
Hii regions in NGC 157 predicted by the models consid- 
ered here, to be compared with the Lyc flux (bottom line) 
required to maintain the DIG in this galaxy ionized, as- 
suming case B recombination and T^IO^ K (Paper I). 



Model / observation 




Lyc phot s ^ 


Density bounding (constant (p) 


14.8 




Density bounding (constant p) 


29.4 




Constant escape fraction 30% 


6.81 




Mixed model 


15.8 


- 30.4 


Required Lyc flux 


12 - 


19 



ized as explained, e.g. in Fig. ^ via histograms showing 
the intensity distribution of all the pixels in these images 
of residuals. This method allowed us to discriminate be- 
tween rather poor solutions (such as the Sciama model 
cited above) or random Ha intensity distributions used as 
trials) and reasonably good solutions, which cover most of 
the cases in which the Lyc photons escape from the Hii 
regions. 

However, it is much more difficult using formal statis- 
tical tests to reach a firmer conclusion about which of the 
assumptions about the escaping fraction of ionizing pho- 
tons yields the best fit to the observations. The models 
which give low residuals are also those that seem to of- 
fer a better result by simple eye estimates, whereas those 
which give higher and more asymmetrically distributed 
residuals are also worse under visual inspection. 

Direct inspection of the figures, and consideration 
of the corresponding residual histograms give sufficient 
grounds to reject the models with the highest values of 
the empirical extinction coefficient: k > 5.3 x 10~^'^ x zioo 
cm~^, while the model which assumes that photons es- 
cape only from the most luminous regions, those with 
Lho > Lstr, does not predict sufficient intensity in the 
DIG near the northernmost arm of the galaxy. Here the 
majority of the regions have Lhol < Lstr but a measurable 
surface brightness is measured in the DIG surrounding the 
arm. It is clear that flux must be escaping from regions 
with Lhci < Lstn a result found also in models where 
extinction was not considered. Very good fits to the ob- 
served DIG surface brightness distribution are found in 
those models where we have taken a constant escape frac- 
tion for the Lyc photons. This appears to go against our 
hypothesis, previously put forward to explain the form of 
the H II region luminosity function in Ha (Beckman et 
al. 2000; Rozas et al. 1996a, 1999, 2000) that Hli regions 
with Lhch > Lstr are the principal sources of ionization 
of the DIG. However our models with a constant escape 
fraction of the ionizing photons down-converted to Ha, 
are normalized and yield maps with only relative intensity 
isophotes in Ha. To fit the total integrated luminosity in 
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the DIG of NGC 157 a constant escape fraction of 70% 
would be required (see Table ^ , and this seems improba- 
ble from a range of physical considerations governing the 
escape of Lyc photons from star forming regions. There- 
fore it seems that if a constant fraction of ionizing photons 
is escaping from the H ii regions, other principal ionizing 
sources (rather than H ii regions) would be needed in or- 
der to explain the power requirements to ionize the DIG 
(Table 0). 

However, fits of very good quality are also found from 
models with a constant escape fraction of ~20% for low 
luminosity regions and an increasing fraction for those re- 
gions with Lhoi > Lstri and this would be in better accord 
with our deductions from luminosity functions and H ii re- 
gion central surface brightnesses (Beckman et al. 2000). 

We should note here, that other ionizing sources may 
contribute to the ionization of the DIG, and the aim of 
this paper is not to demonstrate that the H ii regions are 
the only ionizing source of the DIG. However, the results 
presented give strong support to the hypothesis that H ii 
regions are an important source and probably the main 
source for the ionization. Known possible ionizing sources 
for the DIG in normal galaxies include photoionization 
(by supernovae, Wolf-Rayets, white dwarfs, cosmic rays. 
X-rays, or field OB stars for example, Hoopes & Walter- 
bos 2000), but also shocks and turbulent dissipation in 
the ISM (e.g. Rand 2000; Walterbos 1999) may be im- 
portant in specific zones, as has been noted within the 
Galaxy (Tufte, Reynolds & Haffner 1999). Contributions 
from such sources to the general ionization of the DIG are 
not excluded by the evidence presented here. 

To go further, reliably, into the problem would require 
an effort of detailed modelling of transport in a fully three- 
dimensional clumpy medium, plus the acquisition of an 
H I map of sufficiently high angular resolution. The limi- 
ted conclusion of the present study is that the scenario 
in which leaky H ii regions provide the ionizing sources 
which yield at least the major fraction of the Ha in the 
DIG appears to be well supported in the case of NGC 157. 
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